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Abstract 

Surface waves on liquids act as a dynamical phase grating for incident light. In this article, we 
revisit the classical method of probing such waves (wavelengths of the order of mm) as well as 
inherent properties of liquids and liquid films on liquids, using optical diffraction. A combination 
of simulation and experiment is proposed to trace out the surface wave profiles in various situations 
{eg. for one or more vertical, slightly immersed, electrically driven exciters). Subsequently, the 
surface tension and the spatial damping coefficient (related to viscosity) of a variety of liquids are 
measured carefully in order to gauge the efficiency of measuring liquid properties using this optical 
probe. The final set of results deal with liquid films where dispersion relations, surface and interface 
modes, interfacial tension and related issues are investigated in some detail, both theoretically and 
experimentally. On the whole, our observations and analysis seem to support the claim that this 
simple, low-cost apparatus is capable of providing a wealth of information on liquids and liquid 
surface waves in a non-destructive way. 
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I. INTRODUCTION 



Optical probes are usually non-destructive - they do not change the properties of the 
medium being probed in any significant way (unless the probe power is large). Diffraction 
and scattering of light, as is known over centuries, are therefore capable of providing a variety 
of information. In this article, we revisit a well-known classical method of investigating 
liquid surface waves using light diffraction. Although this method is well-studied in the 
context of measurements of surface tension of liquids [l^ , our motivation for revisiting it is 
primarily aimed at demonstrating its utility and capability in understanding a wider class 
of situations/phenomena. In particular, we investigate two aspects - surface wave profiles 
and liquid properties for waves on a single liquid and also for waves on liquid films. 

Before we begin with details, it is worth summarizing some of the earlier work on liquid 
as well as solid surface waves. Use of laser light scattering is a unique technique for non- 
contact, noninvasive study of surface acoustic waves and capillary waves P|. Technological 
advances directed toward diverse applications of surface acoustic waves on solid surfaces 
are now well established, and has resulted in important developments in the fabrication of 
practical devices 

However, relati vely few articles deal with surface waves on liquid surfaces 
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r, relati vely 
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171 Il8l Il9l . |20[ . To mention a few, high frequency capillary waves at the surface of 
liquid gallium and the mercury liquid-vapor interface have been studied by means of quasi- 
elastic light-scattering spectroscopy 2l| . The spatial damping coefficient of low- frequency 
surface waves at air-water interfaces, using a novel heterodyne light-scattering technique, 
has been measured and reported in [2^. Quasi-elastic light scattering has been used to 
study thermally excited capillary waves on free liquid surfaces over a considerably wider 
range of surface wavenumbers Q| . Measurements of the decay coefficient of capillary 
waves in liquids covered with monolayers of stearic acid, oleyl alcohol, and hemicyanine 
are reported in Also, the measurement of the pressure-area compression isotherm in 
Langmuir monolayer films using the laser light diffraction from surface capillary waves is 
reported in 24 1 . 

In a recent article, we have shown how the effect of interfering waves on a liquid surface 
could be inferred from the nature of the diffraction pattern . Here we begin by proposing 
a combined method of simulation and experiment, sufficiently general in nature and capable 
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of verifying the nature of the actual hquid surface wave profile. Having known the profile, 
we then focus on the dispersion relation and discuss experiments to determine the surface 
tension and viscosity of liquids. Finally, we turn to liquid films where we try to understand 
the profile through the various modes (surface and interface) and then go on to check the 
values of interfacial tension. 

Section II describes the basic scheme of our experimental arrangement. The methodology 
of our experiment and typical results pertaining to the surface wave profile are discussed 
in section III. In section IV, we discuss theory, observations and measurements of surface 
tension and viscosity for different liquids. Section V first discusses the theoretical platform 
for surface capillary wave on liquid films. We then verify our theory with experiments. 
Finally, in section VI, we summarize our results with some concluding remarks. 



II. EXPERIMENTAL DETAILS AND THE THEORETICAL BACKGROUND 



A. Experimental setup 




FIG. 1: Schematic diagram of the experimental setup to study capillary wave on liquid; A: Laser, 
B: Frequency generator, C: Screen, D: Loudspeaker, E: Spectrometer and prism table assembly 
and F: Exciter. 

As shown in Fig. Q a petridish of diameter about 18.5 cm is filled with the experimental 
liquid to a depth of ~ 1 cm. A metal pin with its blunt end glued vertically upright to the 
diaphragm of a loud speaker (held above the petridish) acts as an exciter. When slightly 
immersed in the liquid and driven by a low frequency sinusoidal signal generator, this exciter 
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vibrates vertically up and down and generates the desired liquid surface waves. To study the 
profile of surface waves we place the petridish-loud-speaker assembly on the prism-table of a 
spectrometer in order to choose the diffracting element at any position on the liquid surface 
and to record the corresponding position from the prism-table scale (with a least count of 
0.1 degree). Light from a 5 mW He-Ne laser (A = 632.8 nm) having a beam diameter of 
~ 1.8 mm is directed to fall on the liquid surface at an angle of incidence of 77° {9i), as 
noted in our experiment. The laser beam incident on the dynamical phase grating formed 
on the liquid surface, produces Fraunhofer diffraction pattern which is observed on a screen 
placed at a fairly large distance (3.42 meters in our case) from the diffraction center. The 
images of the diffraction pattern are recorded using a digital camera (make Sony, DSC - 
P93, 5.1 mega pixels). Spurious noise in the image is removed using Photoshop software. 
The intensity across the diffraction spots are measured using a photo-detector operating 
in the linear dynamic range of the photo-currents' response to optical power. During our 
experiment, the room temperature recorded is 25^0. 



B. Theoretical background 



We discuss here, briefly, the diffraction mechanism which lies at the heart of all our 
experiments. When monochromatic light of wavelength A is incident on a circular surface 
wave of frequency u and wavevector K, the phase modulation produced by the surface wave 
is given as 
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cos^j appears due to oblique incidence of the incident monochromatic 
25| |. The field strength E of the diffraction pattern can be estimated from 
the Fourier transform of the aperture function (in this case the surface wave produced by 
sing le exciter). The intensity of the diffraction patterns, obtained from EE*, is given by 

133 

(2) 
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where z is the horizontal distance between the location of the laser spot on the liquid surface 
and the screen and A = is the wavelength of the surface wave, h is the amplitude of the 
surface wave, x' is the coordinate which measures the distance of the diffraction spots from 
a reference point (central spot) on the observation plane. J„ is the Bessel function of order 



n and 5() denotes the Dirac delta function. The intensity distribution on the diffraction 
pattern vanishes at points where J„ is zero satisfying Eqn. 2. The quantity involving the 
Bessel function in this equation gives the magnitude of the intensity of diffraction spots for 
different orders ra, while the positions of the spots are given by the delta function. 
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FIG. 2: Theoretical plot of wave amplitude versus intensity ratio for 0^ = 77 , A = 632.8 nm. 
Solid line, dashed line and dotted line stand for the ratios lo/h, h/Io and lo/h respectively. 

When we generate capillary waves by means of a single pin exciter, the amplitude of the 
wave at the point of oscillation is almost the same as the amplitude of the oscillating pin in 
a low viscous liquid [l^. The intensities of different orders depend only on the amplitude 
of the capillary wave [h), if other parameters like 6i and A are kept constant (Eqn. 2). 
Theoretical plots of wave amplitude versus the intensity ratios, lo/h, h/h and lo/h, using 
Eqn. 2 (keeping 6i and A fixed) are shown in Fig. |21 In our experiment, for a particular 
u, the intensity of different order spots in the diffraction pattern are measured by a photo- 
diode detector when the laser beam is focused along a radial line on the liquid surface. The 
intensity ratios like lo/h and /0//2 are then evaluated. By comparing our experimentally 
measured ratios with those obtained from the theoretical plot (Fig. |21), we estimate the 
average value of h. It is to be noted that, for Jq << /i, the value of this ratio lo/h becomes 
very small and could yield erroneous results. In such cases, we used the ratio Ii/Iq for 
estimation of h (however, we were always below the limit for which this ratio becomes very 
large). Again, by measuring x', we determined the wavelength (A) of the surface capillary 
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waves to be 2.1 mm [for more details see |l5|]. For several exciters placed at different types 
of geometric configurations (eg, along a line, or a polygon etc.), the surface wave patterns 
change, and hence Eqn. 2 must also change. 

III. THE SURFACE WAVE PROFILE 
A. Simulation of the profile 




FIG. 3: The simulated surface capillary wave profiles for (a) single exciter and (b) two exciters. 

Low frequency surface waves are not visible to the naked eye. One may appreciate this 
when we simulate the surface wave profile using realistic dimensions and parameters. In 



a previous study jl5|, we measured the surface- wave parameters (wavelength, amplitude 
at the point of excitation) at a given oscillation frequency (220 Hz). Here, we use similar 
typical values to generate theoretically the distribution of waves on the liquid surface. The 
equations, which we use to simulate the wave profile for the cases of single and double-pin 
exciters are (at a fixed time 't', say t=0) 



MX,y) = hexp (^-S^J{X-x^)' + {Y-y,f^ cos (^\/(X - Xi) V (F - yi)') (3) 
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and 

MX^Y) = hexp (^-S^JiX-x.Y + iY-y.f^ cos (^^^ {X - x^f + {Y - y,f^ + 

h exp (^-5^iX - X2)' + (F - ^2)') cos ( y ^(^ ' ^^f + (Y - ^2)') , (4) 

respectively. Here {xi,yi) and {x2,y2) are the centers of oscillations and 6 is the spatial 
damping coefficient of the liquid. To mimic the intended experiments, we assume that each 
wave has the same frequency u, wavelength A, and amplitude h (at the center of oscillation). 
In our simulations, we use the typical values of A = 2.1 mm, h =1.0 mm and D = 8A mm 
(separation between pin-exciters) 6 is chosen to be 0.235 cm^^ (for water at 220 Hz, 
we obtain this value from experiments discussed later in section IV). The corresponding 
surface-wave proffies estimated for the two cases of interest are shown in Fig. Efa) and 
Fig. ini^b), respectively. Following this recipe one can obtain the surface- wave profile for any 
number of oscillation sources. The simulated proffies enable a better visualization, which, 
in turn, can act as a guideline for their actual experimental verification. 




FIG. 4: The schematic of our experiment, which we have followed to verify the simulated two source 
pattern. Oi and O2 are the two exciters. Solid lines are the hyperbolic loci along the direction of 
destructive interference. 
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B. Experimental results for two source pattern 

The simulated profiles are built on the assumptions made in Eqn. 3 and Eqn. 4. It is 
necessary to verify experimentally whether the assumptions are good enough to represent 
the characteristics of the surface waves. We note that the simulated two source interference 
pattern [shown in Fig. El^b)] exhibits hyperbolic loci for the constructive and destructive 
interference nodes. Using this fact as a guideline, we trace the hyperbolic loci (maxima 
and minima of the interfering waves) on the liquid surface by observing the changes in 
the corresponding diffraction patterns. Initially, when the laser beam is incident along the 
central line PQ (see Fig. pi) where constructive interference of waves occur, the diffracted 
light shows, as expected |l5|, a central spot along with higher order ones symmetrically 
located on either side of the central spot [see Fig. Ell^a)]. As we shift the probe beam 
gradually away from the central line along a circle (note the arrow in Fig. |^, the number of 
higher order spots progressively decreases till we reach the position of the adjacent hyperbola 
(O2A), where the diffraction should ideally contain only the central order. In experiment, 
however, at this position, we find the first order spots too, though with very low intensity 
[see Fig. El^b)]. This happens due to the finite spot size of the laser beam. Nevertheless, 
beginning with the central line at one side and gradually rotating the prism-table till we reach 
180°, we observe a repetition of the patterns [shown in Fig. Efa) andEU^b)] symmetrically 
at all complementary angles. We note the angular positions of all consecutive minima. 
These results are shown in Table-I. Thus, a comparison of our experiment (Fig. |3|) with the 
simulation [upper half of Fig. El^b)] for identical values of parameters shows the same number 
(eight) of destructive interference lines at almost the same angular separation. Beyond 180'' 
and upto 360°, {i.e, while scanning the lower semicircular region in Fig. identical patterns 
are observed. This confirms that the Eqn. 4 used to describe and simulate the two source 
surface wave profile is, indeed, a correct representation of the observed profile. We note that 
our method of investigation is sufficiently general (we have indeed checked upto six exciters) 
and can be used for several exciters located along a line or in other geometric configurations 
(eg. exciters on a triangle, quadrilateral, pentagon, hexagon etc.). 



8 



FIG. 5: Diffraction patterns when light is incident along the direction of (a) constructive interfer- 
ence and (b) destructive interference for the two source (exciter) case. 

C. Minima distribution for oscillation sources in regular array 

When multiple oscillating sources are used, different surface wave profiles can be gener- 
ated by appropriately configuring their positions. It is interesting to study these interference 
patterns of the surface waves for multiple sources (e.g., a linear chain of oscillation sources 
or those placed in a closed regular geometric figure). A typical example of a simulated 
interference pattern for three oscillation sources is shown in Fig. IHl For this purpose, we 
choose to trace the distribution of minima lines in the simulated interference pattern in 
terms of the number of lines and their angular position. It is obvious that the number of 
minima lines increases with the increase in the number of oscillation sources. For example, 
our simulation of the pattern for 2, 3, 4, 5 and 6 oscillation sources placed in a straight 
line with a spacing of D = 8.4 mm between the adjacent two oscillators, and A = 2.1 mm, 
exhibits the number of hyperbolic lines respectively as 8, 16, 24, 32, and 40. To derive 
this behavior theoretically, we use the interference treatment for n-oscillating sources. We 
consider n (> 2) sources lying in a straight line. Representing the wave generated by each 
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oscillating source as hexp{i(K.rn — ujt)), where h is the amplitude and rn is the position 
vector from the n-th source to the point of observation, we find the resultant field ipp a.t & 
certain point P as 

where R = \D{n — l)Cos(0) + ri, is the distance from the centre of the line of oscillators to 
the point P and (3 = KDCos{(f)), is the phase difference between adjacent sources. Again 
is the angle between R and the line joining the sources of oscillation. For minima at P, the 
condition is ^"^"^ly = 0- Using the boundary condition mentioned in [3], we calculate the 
allowed number of minima lines and their corresponding angular positions. For n = 3, for 
example, the above minima condition becomes ^Cos(0) = ±(y), where m is an integer but 
not equal to or multiple of n. If we choose = 4, then Cos(0) has only 16 allowed values, i.e. 
16 hyperbolic minima lines are possible. The angular ditribution of these asymptotic lines 
perfectly match with the observed pattern. We have carried out this comparison successfully 
for higher order n-values upto 5. From the above calculations we have seen that the number 
(No) of asymptotic minima lines in the region above (or below) the line joining the oscillation 
sources is directly proportional to the ratio ^ and also the number of the interval {n — 1) 
between the exciters. Taking care of the ± sign of Cos{(j)), we can express No as 

2D 

No = ^{n-l), (6) 

which gives the same result (for different n nalues) as obtained before. This relationship 
can, just by measuring the number of minima lines for given D and n, estimate A of the 
capillary waves. 

Next we study n (> 3) oscillation sources forming a closed loop shape of an equilateral 
triangle, square, pentagon, and a hexagon. The number {Nc) of minima lines (which are 
lying on ether side of a symmetric axis) in this case can be similarly expressed as, 

2D .n. 

The number of lines for the above stated closed geometrical figures at the same values of 
D and A are 12, 16, 20 and 24 respectively. We have not been able to calculate exactly the 
angular distribution of the minima lines. However, using the simulation mentioned above, 
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□ ipQ 2QD 3nD D I QH 3Dn ADQ 

FIG. 6: Simulated interference profile of surface capillary wave for (a) three exciters in a straight 
line and (b) three exciters at the vertices of an equilateral triangle. 

one can determine the interference pattern due to any number of oscillation sources placed 
along a line or at the vertices of polygons (regular or irregular). Such a study might be useful 
for generating a desired dynamic phase grating aperture for diffraction based micro-photonic 
systems. 

IV. THE DISPERSION RELATION AND LIQUID PROPERTIES USING LIGHT 
DIFFRACTION 

Knowing the surface wave profile is, of course, not enough. A crucial element of any wave 
phenomenon is the dispersion relation, which, obviously, contains within it, quantifiers of 
material properties (here, liquid properties, such as surface tension and viscosity). In this 
section we describe how we use Fraunhofer diffraction of laser light by the surface profile of 
ripples on liquids, described in Section III, to study such liquid properties. Several methods 
to accurately determine these two properties are now well-established and available in the 
literature. For measurement of surface tension, techniques widely used are the Capillary Rise 
method, Drop-weig 

for the measurement of viscosity, popular methods like Poiseuille's method. Stokes' method 
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Incident light 




FIG. 7: Geometrical representation of the basic measurement on sm'face tension and viscosity. 



etc., are well-known and gives accurate results 
set up which can measure both surface tension and viscosity of low-viscous liquids with 
reasonable accuracy. Though the measurement of surface tension of liquids using Fraunhofer 
diffraction of laser light by surface waves is well known in the literature to the best of 
our knowledge the estimation of the viscosity of a liquid using this probe and, in particular, 
the simple method which we have followed, is certainly new. In the following we elaborate 
on the above mentioned measurements in some detail. 



Surface tension 



m 



Consider the well-known dispersion relation for surface capillary wave on liquids |29l . |30 



u^={gK + aKyp). (8) 

Here 00 and K are the angular frequency and wavevector of the surface capillary wave, 
respectively, a and p are the surface tension and density of the liquid and g is the acceleration 
due to gravity. The above relation shows that waves at the surface of the liquid are dominated 
by gravity as well as surface tension. Neglecting the gravitational effect and taking the 
logarithm on both sides of the Eqn. 8 we get 

lnuj = hnK+hn(^^y (9) 

which is an equation of a straight line for Inu vs. InK with a slope of | and y-intercept 
I In (^^^ , which contains the surface tension a. 

In our experiments, we have generated the surface capillary wave on the liquid surface 
with a single exciter. As we discussed before, the expected surface wave profile [shown in 
Fig. Et^a)] acts as a reflection phase grating for the incident laser light and the intensity 
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distribution satisfies Eqn. 2. If we take tlie separation between tlie central order and tlie 
first order (^=1) diffraction spots to be d (see Fig. [7j), tlien from tlie delta function in Eqn. 
2 we get the expression for the wavenumber K as 



K = ^sin^, (10) 
az 

where 6' = ^ — 6*^ is the grazing angle of incidence of the laser beam. For small grazing angle, 
we can approximate K as 



^=T7^' (11) 



2nld 

where / is the perpendicular height of the central order spot from the liquid surface level 
(see Fig. Ej). As changes, the distance between bright spots of the diffraction pattern also 
changes. Sets of diffraction patterns are obtained on the screen for different values of a;. For 
each set, we have traced the positions of spots of different order on the screen. 

For Fraunhofer diffraction, values of z and / in Eqn. 11 are quite large (342 cm and 87 
cm, respectively, in our experiment). Thus, the systematic error introduced in measuring 
these lengths with a metre scale (with a least count of 1.0 mm) is quite less. On the other 
hand, a careful and accurate measurement of d (of the order of few mm) is crucial for 
evaluating K. From the characteristic feature of the delta function in Eqn. 2, we find that 
the diffraction spots due to different orders are equidistant for a particular frequency of 
capillary wave. Thus, to reduce the systematic error in measuring d, we have measured the 
separation between the positive and the negative fourth-order (fifth-order) spots. The values 
of d are obtained through averaging these data. For each frequency u, we have repeated 
the experiment at least four times. Hence, using Eqn. 11, we estimate the value of K for a 
particular u. In addition, it is reasonable to assume that the frequency, of capillary wave 
is the same as that of the exciter (driven by a function generator) [15]. Thus, the values 
of u are found directly from the frequency readings of the function generator driving the 
exciter. Following these measurements, we plot a In -In graph for K vs. uo and evaluate the 
slope and y- intercept of this curve. The slope of the curve (|) should verify the power law 
predicted by the dispersion relation while the y-intercept should give the surface tension. 

Our experimental results on three different liquids (water, 1-hexanol and kerosene - a 
mixture of petrolium hydrocarbons) are shown in Fig. |H| The non-linear least square fit 
(shown by lines in Fig. IH)) to the data with Eqn. 9 verifies the 1.5 exponent of the dispersion 
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FIG. 8: The In-ln plots of the wavevector (K) versus the frequency (w) exhibit a straight hue 
of slope very close to | obtained for three different liquids. The inset shows the corresponding 
y-intercepts. 

relation within experimental uncertainty (1.5±0.03). Using the known value of p for each 
liquid, we estimate the value of a from the measured y-intercept of the straight lines. It is 
to be noted that the maximum percentage error in measuring a is limited to 6% - 8% (for 
different liquids) by our measurement procedure. We have tabulated the measured values 
of surface tension for different liquids in Table-II. These results match quite well with the 
standard values of surface tension for the corresponding liquids j^. 



B. Viscosity 

Using the same experimental set-up, we now move on towards estimating the viscosity 
of a given liquid. We mentioned before that for capillary waves generated by means of a 



lation is almost the same 
However, if we imagine 



single pin exciter, the amplitude of the wave at the point of osci 
as the amplitude of the oscillating pin in a low viscous liquid ^ 
a radial straight line from the source of excitation (on the liquid surface), then along this 
line the amplitude of surface wave is gradually damped due to the viscosity of the liquid 
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FIG. 9: Experimental plot of distance (x) versus surface capillary wave amplitude (h). Solid line, 
dashed line, and dotted line for the liquids water, kerosene, and 1-hexanol respectively. 

(see Fig. The laser beam is focused along a radial line at different distances on the 
liquid surface from the point of excitation {x in Fig. [7j). We have seen in Sec. II that the 
intensities of the diffraction spots of different orders depend only on the amplitude of the 
capillary wave, if other parameters like 6i and A are kept constant (Eqn. 2). By comparing 
our experimentally measured intensity ratios of different orders in the diffraction pattern 
with those obtained from the theoretical plot (Fig. we estimate the average value of h 
at a particular location on the liquid surface. In Fig. IHl we plot x, the distance from the 
point of oscillation vs. h, the capillary wave amplitude for three experimental liquids (water, 
kerosene, and 1-hexanol). We fit the data points with an exponential decay function, where 
we kept the decay constant (also known as spatial damping coefficient of the liquid), 6, as 
varying parameter. The values of 6 for three different liquids are tabulated in Table-II. 

The wavevector for a wave of fixed real frequency, u, has a small imaginary part, which 
contributes to the spatial damping. Substituting K = KQ + i6 (where 6 « K g) in the well- 



known Navier-Stokes equation for waves on a liquid-air interface, one gets|22l. 



33| 



(12) 



3a ' 

where t] is the viscosity of the liquid. Experimentally, we have obtained the values of a and 
6 for the above mentioned liquids. The estimated value of viscosity of different liquids, using 
Eqn. 12, are shown in Table-II. The experimental plots in Fig. El are for frequency 220 Hz. 
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We have also checked our results for two other frequencies (260 Hz and 300 Hz) and have 
obtained similar results (see Table-II). 



V. LIQUID FILMS ON LIQUIDS: THEORY AND EXPERIMENTS 

Till now we have been exclusively concerned with waves on the surface of a single liquid, 
or, more precisely, waves at the liquid-air interface. A more complicated scenario arises 
when there exists a film of a different liquid on top of a given liquid (the two liquids are 
immiscible). In such a situation, the dispersion relation changes drastically, which, in turn, 
changes the surface wave profile. We now investigate the novelties that arise through studies 
with our optical probe. 



A. Background theoretical framework 

The dispersion relation in Eqn. 8 for capillary wave propagation on liquid surfaces is 



valid only when the depth of experimental liquid is 
the liquid (denoted by, say y), the relation becomes 



arge. For lower values of the depth of 

m 



uj^ = ^tanh{Ky) (13) 

where, here too, we neglect a term due to gravity. Typically, if we choose y = 1 cm and a 
range of K from 20 to 80 cm ~^ (as used in our earlier experiment), tanh {Ky) is nearly unity, 
and the dispersion relation reduces to the one given by Eqn. 8. The interesting problem we 
focus on now is that of wave propagation on a liquid film spread out on another liquid. If the 
thickness of the film is very low, then not only the surface tension but also the interfacial 
tension across the liquid-film boundary, play a role in the propagation mechanism of the 
waves. We first work out a general dispersion relation where the surface and interfacial 
tensions as well as the film thickness appear explicitly in our mathematical expressions. 

Fig. ITUl shows the symbols used to represent various parameters in the following discus- 
sion. We denote the density of the lower liquid by p while a is the interfacial tension at the 
liquid-liquid interface. The same quantities for the film- air interface are represented with 
primes i.e., p and a respectively. The equilibrium plane of separation (interface) between 
the liquids is at ?/ = 0. The liquid below the y = plane extends upto a value y = —H 
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FIG. 10: Liquid film on liquid and the various parameters used. 

with H being very large. The film is located above y = with a thickness y' << H. Under 
these assumptions, the general expressions for the surface capillary waves in the region (I): 
-H < y < and in region II : < y < y' are given by 



^ = Ae^y cos {Kx - cot) 



(14) 



and 



Be-^y + Ce^y\ cos {Kx - uot) , (15) 

respectively. Here A, B and C are arbitrary constants. The above two waves must obey 
the following boundary conditions arising from the continuity of the velocity component (fy) 
and that of pressure at the surface of separation of the two liquids, i.e., 

dip dip' 



a 



P 



'd^iP' 



P 



— a 



, d ( d'^ip' 



9P 



0. 



'dip' 
, dy 



(16) 
(17) 

dy ' ^ dt"^ dy \ dx^ J ^ 
Evidently, the first two conditions apply at y = while the third one is valid at the top 

surface aX y = y . In the boundary conditions above, we neglect the effect of gravity in all 

calculations below. Using Eqn. 14 and Eqn. 15 in Eqn. 16 and Eqn. 17 we obtain 



dtp d'^ip 
''Ty + 

, dip' 
P9- 



dy dy 

d_ (d^ 

dy \ dx^ ^ 
,d^iP' 



and 



C = B 



A = C -B 



aK^ - (p + p') 



[aK^ — {p — p')uj'^ 



(19) 



(20) 
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Substituting Eqn. 19 and Eqn. 20 in Eqn. 18 with y = y' ,we construct a quadratic equation 
for u;^ given by: 

[l + r tanh (^2/')] - oo^ [^^ tanh (fsT?/') + rw^ + cjf] + cu^cuf tanh (fsT?/') = 0, (21) 

where, we have used the notations ^ = r, = u"^ and = ujf for simphcity. The two 
roots of Eqn. 21 are the new dispersion relations for our given problem. The two roots are 

4 = (22) 

where, P = ruj^ + ujlia.-nh.{Ky'^ + oof, Q = 1 + rtaiih.(^Ky''j, and G = P'^ — 

Consider two limiting cases of the above dispersion relation. Let us first assume y ^ oo. 
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Then u"]^ —>■ oOg and cj^ = Thus, when y' is very large, the mode corresponds 

to the waves propagating on the upper liquid surface (thus called surface mode). On the 
other hand, the mode involves the interfacial tension as well as the density of the 
lower liquid and may be termed as an interface mode. In such a situation {y' large), reflecting 
the probe laser beam off the surface wave profile on the liquid film, we can obtain information 
(through the diffraction pattern) on wave propagation and properties of the liquid film. The 
existence of the other mode can be confirmed (for y' large) only if we perform experiments 
in transmitted light, though it is not clear how such information may be obtained. On 

the contrary, assuming ?/' — ^ (very thin film), we get cj^ — > [rujl + ujf) = - — — , and 
u"^ — * 0. Thus the mode now depends on properties of both the liquids, {a', p), as well 
as properties of the interface (a). 

Furthermore, we can obtain the constants A, B, C and hence i}) and if)' by assuming a 
value for the amplitude imparted initially at the top surface. In this way, one may obtain 
the amplitudes of il) and ip' for both the modes uo^ and uj^. In particular, we have noted (not 
demonstrated here) that the nature of variation of the amplitudes of ip' (with the film-hight 
< y' < h) for the u;+ mode is quite different from that for the uj_ mode. In the former case 
iuj^) the value of the amplitude drops to a much smaller value at the interface compared to 
the drop in the amplitude for the u;_ mode. However, it is not possible for us to check this 
fact through our experiments. Hence we refrain from discussing this aspect further in this 
article. 
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B. Experiments and results 




FIG. 11: Theoretical plot of frequency (w) versus wavevector (K) for different modes in water-o- 
xylene liquid- liquid film at film thickness (a) 0.158 cm, (b) 0.016 cm and (c) 0.032 cm. The solid 
line is for u+ mode, dashed line is for co- mode and dotted line is for cOg mode. The experimental 
data points are shown by the filled circles. 

In our experiments, we have tried to figure out which of the two different modes (discussed 
in previous section) that arise due to the presence of films of kerosene, n,n-dimethyl aniline, or 
o-xylene spread on the surface of water, is actually present. The choice of these liquids is 
due to the fact that these are immiscible with water, have density lower w.r.t. water, are 
low- viscous and have a higher boiling point. In these experiments with films, a petridish of 
area ~ 63 cm^ is filled with water to a depth of nearly 1 cm and the liquid film is formed 
by adding a measured quantity of the second liquid in drops from a micropipctte. The 
thickness of the film is estimated from the known volume of the liquid and the surface area 
of the petridish. Uniformity of the films on the water surface has been checked optically. 
In order to avoid any decrease in thickness due to evaporation, we have taken observations 
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within 5 minutes of the formation of the uniformly spread film. The measurements are of 
the same type as those mentioned while studying surface tension (discussed in the earlier 
section) except that, in the present case, we do the experiments for varying film thicknesses. 
Here, we show our experimental results only for o- xylene films of different thicknesses on 
water. For films of the other liquids mentioned above, on water, we have observed similar 
behavior. 

To find out the mode which is present on the film surface for a liquid-film of o-xylene on 
water, we first plot the experimental uj ^ K behavior with film thickness as a parameter. 
The theoretically obtained uo-K curves for different modes : cu-, and uOg (as a limiting 
case of ujj^ mode) - are then placed on the same graph for comparison. Here we use standard 
values : a = 37.2 dyne/cm, a' = 29.8 dyne/cm, p'= 0.88 gm/cc, and p=l gm/cc in Eqn. 

cm, the measured data 

corresponds to mode [Fig. [TlT b)]. while for films of thickness typically 0.158 cm, the 
experimental data matches well with the Ug mode [Fig. fTlT a)]. However, any intermediate 
thickness between 0.016 and 0.158 cm yields data within the u+ and ujg modes [Fig. ITTT c)]. 
This clearly shows that the Uj^ mode always dominates over U- mode. For thicker films the 
0;+ mode coincides with Ug mode, whereas for thinner films the mode is still though it 
is modified by the presence of the interfacial tension and the density of the lower liquid and 
is not equal to Ug. 

The above conclusions have been further verified by analyzing the experimental data for 
the estimation of the surface and interfacial tension. We have fitted our experimental data 
of water-o-xylene film of thickness 0.016 cm with 00+ and uj- modes in Eqn. 22 by nonlinear 
curve fitting (keeping surface tension and interfacial tension as varying parameters). The 
UJ+ mode fits with the experimental data and estimates the values of surface tension [a') 
of o-xylene as (29.5 ± 3.6) dyne/cm and the values of interfacial tension (a) as (33.2 ± 8.8) 
dy„e/c,„ Isolid linein Fig. ig. These values are very dose to the standard values available 
in the literature [2^. On the other hand, uj- mode fails to fit the experimental data well 
(dashed line in Fig. fT^ : moreover, it yields unphysical values of surface and interfacial 
tension. Experimental limitations, in our present set up, did not allow us to make films of 
thickness lower than the values quoted above. 
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FIG. 12: Nonlinear curve fitting of the experimental data for water- o-xylene film of thickness 
0.016 cm with a;+ and U- modes. 

VI. CONCLUDING REMARKS 

In this article, we describe a simple optics-based experimental technique devised with the 
necessary background theoretical formulation to study the characteristics of surface capillary 
waves on liquids. Firstly, we visualize the profile of the surface waves using simulations with 
realistic values for the parameters involved. The locations of the nodes and antinodes which 
appear in the simulated profile are used to match the profile with actual experiments. We 
also claim that this combined method of simulation and experimentation can be used for 
other geometric configurations of oscillating sources on liquid surfaces. Subsequently, the 
two key liquid properties, surface tension and viscosity, which appear in the dispersion 
relation and play an integral role in the propagation of such capillary waves on a liquid 
surface, are measured optically using Fraunhofer diffraction of laser light. We note, in 
particular, the novelty and simplicity in the measurement of viscosity and the fact that 
our measurement tools, are non-destructive but by no means, less precise. Results of our 
experiments with several liquids, when compared to the respective values known from other 
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sources and literature, summarily establish the efficacy and accuracy of our approach. We 
anticipate, in future, the introduction of technological sophistications in our set-up, which 
might lead to an 'Optical Surfacetensometer cum Viscometer'. The advantage of having a 
single set-up for measuring both these properties need not be further emphasised. Finally, we 
investigate, both experimentally and theoretically, an interesting aspect of capillary wave 
characteristics involved with a liquid film placed on the surface of an immissible liquid. 
The dispersion relation we frame to study this case not only explains the features of the 
surface and interfacial modes, but also provides an estimate of the interfacial tension across 
the film-liquid boundary. In summary, our experimental and theoretical results seem to 
demonstrate the capability of our set-up in carrying out studies pertaining to capillary wave 
profile, liquid properties and waves on liquid films on liquids. Despite its limitations, the 
wide variety of information which can be obtained from such a simple set-up, makes it a 
tool worth improving upon in future investigations. 
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Descriptions 


Angle between the lines 




A - B 


B - C 


C - D 


D - E 


E - F 


F - G 


G - H 


Simulated angular differences 
between the two corresponding 


25.5° 


18° 


16° 


15° 


16° 


18° 


25.5° 


minima 
















Experimental angular differences 
between the two corresponding 


25.1° 


17.6° 


16° 


15° 


15.7° 


17.5° 


25.3° 


minima 

















TABLE 1: The simulated and experimental angular separations between the two successive de- 
structive minima of two source interference pattern [see Fig. 2(b) and Fig. 0]. 



Liquids 


a in dyne/cm 
(experimental) 


a in dyne/cm 
(standard) 


b in cgs 
units 


T] in cp 
(experimental) 


T] in cp 
(standard) 


Water 


75 ± 5.3 


72 


(0.235 ± 0.026) 
(0.290 ±0.024) 
(0.304 ±0.031) 


(0.92 ±0.1) at 220 Hz 
(0.96 ± 0.08) at 260 Hz 
(0.87 ± 0.09) at 300 Hz 


0.89 


Kerosene 


26.3 ± 3.7 


28 


(0.947 ±0.04) 
(1.11 ±0.05) 


(1.44 ± 0.06) at 220 Hz 
(1.43 ± 0.07) at 260 Hz 


1.40 


1-Hexanol 


28.5 ± 3.4 


26 


(2.75 ±0.17) 


(4.27 ± 0.26) at 220 Hz 


4.58 



TABLE II: The experimental and standard values of surface tension {a) and viscosity (V) of different 
liquids at room temperature (25° C). For standard values of a and r/ see reference 32| . 
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